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Abstract 

The Clar number of a (hydro)carbon molecule, introduced by Clar [E. Clar, The aromatic sextet, 
(1972).], is the maximum number of mutually disjoint resonant hexagons in the molecule. Calculating 
the Clar number can be formulated as an optimization problem on 2-connected planar graphs. Namely, 
it is the maximum number of mutually disjoint even faces a perfect matching can simultaneously alter¬ 
nate on. It was proved by Abeledo and Atkinson [H. G. Abeledo and G. W. Atkinson, Unimodularity 
of the clar number problem, Linear algebra and its applications 420 (2007), no. 2, 441-448] that the 
Clar number can be computed in polynomial time if the plane graph has even faces only. We prove 
that calculating the Clar number in general 2-connected plane graphs is NP-hard. We also prove 
NP-hardness of the maximum independent set problem for 2-connected plane graphs with odd faces 
only, which may be of independent interest. Finally, we give an FPT algorithm that determines the 
Clar number of a given 2-connected plane graph. The parameter of the algorithm is the length of the 
shortest odd join in the planar dual graph. For fullerenes this is not yet a polynomial algorithm, but 
for certain carbon nanotubes it gives an efficient algorithm. 
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1 Introduction 

Our research is motivated by problems in chemical graph theory. Some molecules, for example polycyclic 
aromatic hydrocarbon (PAH) molecules, benzenoid hydrocarbon molecules, or fullerene molecules can be 
represented as a 2-connected plane graph. In this representation only carbon atoms are depicted, while 
hydrogen atoms are omitted. Several of the chemical properties of these molecules (e.g. chemical stability) 
are closely related to the parameters of the underlying graph. In this note we will be concerned with one 
of these parameters, the Clar number that we will define later. 

A subclass of PAHs, the benzenoid PAHs have the special property that every face in this plane graph 
is a hexagon, in particular, every face has an even number of nodes. One can see that such a 2-connected 
plane graph is also bipartite. Abeledo and Atkinson proved that the Clar number of a 2-connected bipartite 
plane graph can be computed in polynomial time. On the other hand we show in this note that determining 
the Clar number of a general 2-connected plane graph is NP-hard. 

Fullerenes are carbon molecules with a hollow cage-like structure. The first fullerene molecule to 
be discovered, and the family’s namesake, buckminsterfullerene (Ceo), was prepared in 1985 by Richard 
Smalley, Robert Curl, James Heath, Sean O’Brien, and Harold Kroto at Rice University [5]. The graph 
representing a fullerene molecule contains exactly 12 pentagon faces, and the rest of the faces are hexagons 

*Department of Operations Research, Eotvos University, Pazmany Peter setany 1/C, Budapest, Hungary, H-1117. Sup¬ 
ported by the Hungarian Scientific Research Fund (OTKA, grant number K109240). 

^MTA-ELTE Egervary Research Group, Department of Operations Research, Eotvos University, Pazmany Peter setany 
1/C, Budapest, Hungary, H-1117. Supported by the Hungarian Scientific Research Fund (OTKA, grant number K109240). 


1 



(the number of hexagons can be arbitrarily large). For the Clar number of fullerenes Ye and Zhang gave 
an upper bound of [9]. Later they characterized the fullerenes achieving this bound [8]. 

One of our motivations was to determine the Clar number of fullerene molecules in polynomial time. 
We show that determining the Clar number of a general 2-connected plane graph is NP-hard, if the number 
of odd faces is not bounded in the planar embedding We present an algorithm that determines the Clar 
number of a 2-connected plane graph, and has good running time, provided that the odd faces are “not 
too far from each other”. More precisely, our algorithm is fixed parameter tractable (FPT) where the 
parameter is the length of the shortest odd join in the planar dual graph. In Sectionwe explain that for 
a subclass of fullerenes (called carbon nanotubes) our algorithm efficiently computes the Clar number. 

Let G = (V, E) denote a 2-connected planar graph which has a perfect matching. For a planar 
embedding of G and a perfect matching of G let Fm denote the set of those faces which alternate with 
respect to M. Note that faces in Fm are even. A pairwise vertex disjoint subset of Fm is a Clar set with 
respect to M. A subset G of the faces is a Clar set if there exists a perfect matching M for which G is 
a Clar set with respect to M. Note that a set of pairwise vertex disjoint even faces is a Clar set if and only 
if deleting all (the nodes of) these even faces the remaining graph still has a perfect matching. The Clar 
number of G, denoted by Gl{G) is the maximum size of a Clar set. For sake of simplicity we allow the 
unbounded face in a Clar set as well, but there are no difficulties if we want to exclude it. By a plane 
graph we mean a planar graph with a fixed planar embedding. For further graph theoretic definitions we 
refer the reader to [4j. 

The Clar number was defined by Clar in [3]. It was proved by Abeledo and Atkinson [2] that the Clar 
number can be computed in polynomial time if G is bipartite. Note that a plane graph is bipartite if and 
only if all its faces are even. Our first result is the following theorem. 

Theorem 1. It is NP-hard to calculate the Clar number of a 2-connected planar graph (given with a fixed 
planar embedding). 

The proof of the above theorem will be detailed in Section based on a reduction from a special case of 
the independent set problem. In Section]^ we present an FPT algorithm that determines the Clar number 
of a 2-connected plane graph. 

2 Hardness of the Clar number problem 

In this section we prove Theorem Our reduction will be based on a special case of the Independent Set 
Problem. Let us start with defining this problem. 

Definition 1. Given a graph G = (V,E), a subset U G V is said to be independent if there is no edge 
of G between two nodes of U. Let a{G) denote the maximum size of an independent set in G. 

Problem 1. Given a 2-connected planar cubic graph G and a positive integer K, does G contain an 
independent set of size K ? 

Theorem 2 (Mohar, Theorem 4.1 in [5]). Prohlem^is NP-complete. 

Problem 2. Given a 2-connected plane graph G with odd faces only, and a positive integer K, does G 
contain an independent set of size K ? 

Lemma 1. Problem\^is NP-hard. 

Proof. According to Theorem the independent set problem is also NP-complete for 2-connected planar 
graphs. Let G = {V, E) denote an instance of this problem, and let us fix a planar embedding of G. If G 
has an even face F, let Gp denote the planar graph obtained from G by the following operation. We add 
three vertices a, b, c inside F and edges ab, be, ca, au, bu, bv where u and v form an edge of F (see Figure 

0 ^ _ 

^Fullerene graphs are quite special 2-connected plane graphs: for example they are 3-connected and 3-regular, too. We 
believe that among these properties the constant bound 12 on the number of odd faces is the most important, if we are 
concerned with their Clar number. 
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Figure 1: Eliminating even faces. 


Claim 1. a^Gp) = 1 + a{G). 

Proof. First, for an independent set I of G, clearly lUc is independent in Gp and hence a{Gp) > H-a(G). 
Second, an independent set Ip in Gp can contain at most one vertex from the set {a, 5, c}. Since Ip\{a, 6, c} 
is independent in G we get that a(G) > a(Gp) — 1. □ 

Note that the number of even faces of Gf is one less than that of G, and Gf is also 2-connected. Let 
F denote the set of even faces of G. By consecutively applying the above operation on every member of F 
we get another graph Gp for which a(Gr) = a(G) + |F| and which has odd faces only. □ 

We are now ready to prove the hardness of the Clar number problem. 

Proof of Theorem^ We prove the theorem by reducing Problem to the Clar number problem. Let 
G = {V, E) denote an instance of this problem. We construct graph G' the following way: for every edge 
of G we add two vertices to G'. Let uv € E he an edge of G and let Ex and denote the faces uv is 
incident to. We add vertices Xuv,Fi and Xuv,F 2 to G' along with the edge Xuv.FiXuv.f^- If edges uv and vw 
are neighbouring edges on a face F, then we add edge Xuv,fXvw,f to G'. It is easy to see that G' is planar 
(see Figure]^. Informally, G' is obtained from the planar dual graph G* of G by “blowing a circuit” into 
each vertex of G*. Every face of G' either corresponds to a face of G, or a vertex of G, and since G has 
odd faces only, all the even faces of G' are the ones corresponding to vertices of G. Note that G' trivially 
has a perfect matching M consisting of the edges of the form Xuv,FiXuv,F 2 ^ for every uv G E. Since M is 
alternating on every even face of G', corresponding to a vertex of G, that is, on every even face of G', for 
this graph the Clar number equals the maximum size of a Clar set with respect to M. The Clar sets of G' 
and the independent sets of G have a one to one correspondence, proving the theorem. □ 

Corollary 1. It is also NP-hard to find a maximum cardinality Clar set with respect to a fixed perfect 
matching. 


3 An FPT algorithm for determining the Clar number 

In this section we present an algorithm that determines the Clar number of a 2-connected plane graph, 
and has a good running time, unless the odd faces are “far from each other” in the planar representation. 
The idea is the following. Consider a 2-connected plane graph that has only 2 odd faces in its (fixed) 
planar representation, and take a shortest path (in the planar dual graph) between these odd faces. An 
optimal Clar set might use some of the even faces that lie on this shortest path. Our algorithm takes an 
arbitrary subset of even faces along this shortest path and tries to extend this subset into a Clar set. This 
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Figure 2: Reduction of the Independent Set Problem to the Clar number problem 


is repeated for every possible subset of even faces along the shortest path. We will generalize this for plane 
graphs having more than 2 odd faces below. First we need a definition and a theorem. 

Definition 2. Given a graph G = {V, E) and a subset T C V of even size, a T-join is a subset of edges 
J f- E so that the number of edges of J incident to a node v G V is odd if and only if v GT . An odd-join 
of G is a T-join where T is the set of nodes having odd degree in G. 

Theorem 3 (See e.g. [7], Chapter 29 ). Given a graph G = {V,E), a subset T C V of even size, and 
edge-lengths c : E ^ K_|_, a shortest T-join can be found in polynomial time. 

Given a 2-connected plane graph G = {V,E), let G* = {V*,E*) denote its planar dual. Let T C V* 

be the set of odd-degree nodes of G*. Let J* C E* he a. shortest T-join in G*, where each edge of G* 

has length 1 (we refer the reader to [3 Chapter 29] for preliminaries on T-joins). We give an algorithm 
determining the Clar number of G that runs in 0(31“^ ^p(l^l)) for some polynomial p. 

Let J C E he the set of edges corresponding to J*. Let Te„en be the set of even faces of G (that 

is, |Te„e„| = |y*| — |T|) and let Ej C Ff,yen be the set of even faces that have some edge of J in their 
boundary. Let G' = {V -\- U, E — J -\- J') be the 2-connected bipartite plane graph that is obtained from 
G by subdividing each edge of J with a new node (the set of these subdivision nodes is U, the set of 
subdivided edges is J': observe that \U\ = | J| and | J'j = 2| J|). Note that G' is indeed bipartite, since 
every face is even in its planar embedding. 

Let K' be the node-edge incidence matrix of G', and R' be the node-face incidence matrix of G'. Let 
R be obtained from R' by deleting the columns corresponding to odd faces of G, and let K be obtained 
from K' by deleting the columns corresponding to J'. Since the matrix [R', K'] is unimodular by Theorem 
3.5 of [T] we get the following claim. 

Claim 2. The matrix [R,K\ is unimodular. 

After these preliminaries we present the pseudocode of our algorithm that calculates the Clar number 
of a 2-connected plane graph. The basic idea of the algorithm is the following. Determining the Clar 
number of G means that we want to choose pairwise node disjoint even faces and edges, so that every node 
is contained in exactly one of the chosen objects, and we want to maximize the number of faces chosen. 
Given a feasible solution consisting of a set Fi of even faces of G and a set Ei of edges of G, let F{ be the 
set of even faces of G' corresponding to faces in Fi , and similarly let E[ be set of edges of G' corresponding 
to edges in Ei (if an edge e € J is in Ti then we add both edges obtained from the subdivision of e into 
E[). Every subdivision node (that is, node in U) is then incident to either 0,1 or 2 objects in F{ U E[. 
If someone tells us these 0,1,2 values for every u G U then we can reconstruct Ei and Fi using these 
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numbers, see Lemma below. Therefore what we do is that we try every possible vector in {0,1, 2}^ to 
find the one giving the best solution. 

Algorithm Clar_Number 
begin 

INPUT: a 2-connected plane graph G= {V,E) 

OUTPUT: the Clar number of G 

1.1. Find a shortest odd-join in G* , where each edge of G* has length 1 {G* is the planar dual of G, and we 
will use more notations that were introduced above in this section). 

1.2. For every vector bjj G {0,1, 2}^ 

1.3. Let ^ e {0,1,2}'"+'^. 

1.4. For every e G J' 

1.5. Let Ze = 1 if e is incident with a node u gU with bu{u) = 2, and let Zg = 0 otherwise. 

1.6. Take the integer optimum of the LP Problem Q-([^ (see Lemma [^. 

maxjly : y G G , (1) 

Ry + K'x = b,Xe = Zg for every e G J'.} (2) 

1.7. Output the best of the candidates obtained in Step |1.6| 
end 

Lemma 2. The LP Problem has an integer optimum. 

Proof. After eliminating the variables Xg for e G J' we obtain an LP Problem of the form maxjly : y G 
a; G Ry + Kx = 6'}. The polyhedron in this problem is integral by Claimj^ □ 

Note that the LP Problem 0-® will not necessarily be feasible for every choice of bjj. We could 

be more careful in choosing only those vectors in Step |1.2| of the algorithm that make the LP feasible. 
However the algorithm is easier described this way. The running time is clearly 0(3^'^ for some 

polynomial p. 

Carbon nanotubes are fullerenes with a cylindrical nanostructure, with two ’half-fullerene’ caps on both 
ends. Six pentagonal faces are in both caps, forming three short pairs in the odd join. So for this class of 
fullerenes the parameter of our FPT algorithm is relatively small, giving an efficient method to determine 
the Clar number. 


4 Open questions 

We have proved the NP-hardness of the Clar number problem for general plane graph G. The problem 
is motivated by the problem of determining the Clar number of fullerene graphs, when G has exactly 
twelve pentagonal faces and every other face is a hexagon. This problem is however left open, since our 
NP-hardness reduction involves creating a lot of odd faces. An FPT algorithm with the number of odd 
faces as parameter would yield a polynomial time algorithm for all fullerenes. 

Another line of research would be to show that determining the Clar number is NP-hard even for some 
restricted class of 2-connected graphs, too. If we were able to specialize the Independent Set problem 
further to 3-regular plane graphs with odd faces, then our techniques would yield that the Clar number is 
NP-hard for graphs with only hexagonal even faces. 
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